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ABSTRACT 


The charge distribution of periodic electron beam pulses generates patterns 
of sub—Cerenkov radiation distinctive of the distribution of charge with in a 
bunch. Mapping the radiation pattern from different charge shapes may provide 
insight into whether charge pulse shapes can be determined from observed 
radiation patterns. The radiation patterns of Gaussian, Level, and Trapezoidal 
function were mapped by computer simulation. Near 90° to the beam, the 
radiation patterns of all three charge distributions developed an envelope 
proportional to the fourier transform of the charge bunch distribution when the 
wavelength of the emitted radiation was comparable to the size of the bunch. 
For the Gaussian function, the envelope is Gaussian, for the level function it 
is a sinc function. Since the envelope for the trapezoidal function in the 
product of two sinc functions it is more difficult to analyze. This work may 
provide a basis for determining the charge shape of electron beam pulses from 


the sub-Cerenkov radiation based on the radiation intensity pattern. 
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I. INTRODUCTION 


A. HISTORY 

A very faint emission of bluish-white light, what was later-on called 
Cerenkov radiation, from transparent substances placed nearby a radioactive 
sources had been observed by Mme Curie in 1910 . Although electromagnetic 
theory had been sufficiently developed at that time to _ describe this 
phenomenon, the first deliberate attempt to study the phenomenon was made 
by Mallet after many years, who found that the light emitted from a wide 
variety of transparent bodies placed close to radioactive source always had the 
same bluish-white light quality, and that spectrum was continuous, not 
possessing the line or band structure characteristic of fluorescence [Ref. 1: p. 
1]. 

P.A. Cerenkov carried out a series of experiments between 1934 and 1938 
related to the phenomenon. In 1937 [Ilya Frank and Igor Tamm proposed a 
satisfactory theory of the radiation. The experimental results and theoretical 
predictions were in excellent agreement. A complete description of Cerenkov 


radiation, as well as exact mathematical treatment is given in [Ref. 1: p. 1]. 


B. WHAT IS THE CERENKOV RADIATION ? 

A charged particle or a group of charged particles moving through a 
medium at a speed greater than the (phase) velocity of light in the same 
medium will generate a continuous spectrum of electromagnetic radiation. This 


process is referred to as Cerenkov radiation.[Ref. 2: p.9] 


C. PREVIOUS RESEARCH AT N.P.S. 

The research on Cerenkov radiation has been conducted by Professors Fred 
R. Buskirk, John R. Neighbours and others in the Naval Postgraduate school 
since the LINAC was established. Recent work has concentrated on the 
determination of electron beam bunch profile through Cerenkov radiation. 

The work was aimed to illustrate examples of the consequences of the 
finite spatial charge distribution to show that Cerenkov radiation emitted at 
low microwave frequencies may be used to characterize the properties of an 
intense relativiv*; beam in air. The results were that the spatial charge 
distribution of an electron pulse, along with the beam interaction length, 
determines the Cerenkov radiation distribution as a function of frequency. An 
angular distribution of the Cerenkov radiation can, in principle, measure its 


spatial charge distribution.[Ref. 3:p. 1994]. 


D. PURPOSE 

The charge distribution of a bunch influences the pattern of radiation 
produced. Detailed calculations for several simple charge distributions have 
been carried out for the Cerenkov case [Ref. 4:pp. 1994,1996]. The object of 
this work is to assess the effect of the charge distribution of a bunch on the 
sub-Cerenkov radiation from a charged particle beam. Specifically, this means 


investigation of the short beam path region at energies below the threshold. 
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Il. THEORY 


A. BACKGROUND 

A fast, charged particle may cause electromagnetic radiation in a dielectric 
medium when it moves with a constant speed. Suppose an electron to be 
moving relatively slowly through a transparent medium as shown on Fig.1.1 
a). As the particle passes through the medium, each elemental region of the 
medium along the track will in turn receive a very brief electromagnetic pulse 
and be polarized, producing a dipole. An electromagnetic pulse is generated by 
the formation and disappearance of the dipole as the bunch passes. The 
electromagnetic pulse propagates away from the dipole source at the velocity of 
light in the medium. The distance between wave fronts is compressed as the 
velocity of the electron bunch increases. This effect is similar to the doppler 
effect for sound radiated from a moving source. If the speed of the particle is 
less than the propagation velocity in the medium, the electromagnetic field 
produced is called sub—Cerenkov radiation. [Ref. 1: pp. 3-4] 

If the electron is moving faster than light in the medium (see Fig 1.1 b), 
the the pulse source velocity is greater than that of the propagating 
electromagnetic wave which transport the energy. The result is a wavefront 
and where the emanated wavelets bunch together in phase. This shockfront of 


electromagnetic radiation is Cerenkov radiation. [Ref. 1: pp. 4-5] 


1] 





Fig. 2.1 Polarized Atoms in a Dielectric. 


B. CHARACTERISTICS OF CERENKOV RADIATION 
Fig. 2.2 shows the resultant Cerenkov radiation plane wavefront BC from 
particles traveling from A to B. Wavelets from point 0 to 5 are coherent 


along wavefront BC. The angle @ is given by the "Cerenkov relation": 


Ce ome — ee 6) 0 (Qeiy 


where n is the refractive index of the medium and 3 is the ratio of the 
velocity of the particle to the velocity of light in a vacuum [Ref. 1: pp.4—5]. 
From the above, the following is seen that: 
(1) For medium of given refractive index n, there is a threshold velocity 
bac = (1/n), below which no radiation take place [Ref. 1: p. 5]. 
(2) For ultra-relativistic particles, for which 3 = Il.there is a maximum 
angle of emission, given by Omax = cos '(1/n). 


(3) Radiation occurs mainly in the visible and near visible region for 


which n > 1. In the A-ray region, above atomic resonances, i Is 


usually less than 1; therefore Cerenkov radiation in the X-ray region 


is not probable [Ref. 3: p. 3252, Ref. 1: p. 6]. 











-_ 
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ig. 2.2 Radiation Shockfront 


(4) In a gas, @ has only a slight dependence on electron velocity because 
3 must be close to one in order to obtain the Cerenkov condition 
(Ref. 3: p. 3249]. 


From the nature of the radiation, which has been described above, it may be 


I3 


concluded: 
o, Azimuthal symmetry of the radiation produces a cone of radiation. 
(6) The distribution in @ of the radiation intensity approximate a Dirac 
é-function centered about Cerenkov angle ' 0’. 
(7) Length AB in Fig.2 must be much larger than the radiated 
wavelength in order to avoid diffraction effects. For a finite radiation 


length, the diffraction effect is significant. 


C. CALCULATED POWER OF COHERENT CERENKOV RADIATION 
The total coherent power per unit solid angle, radiated at the frequency v 
by a periodic charged particle beam in traveling a finite distance L at 


constant velocity is 


Wik) = ap Oar. (29) 


where vo is the fundamental frequency of the beam generator and v is a 
harmonic of vp, q is the charge of an individual bunch, Q is a constant and R 
is the radiation function. The bunches in the beam are assumed to be 
unchanging in shape and size as the beam travels through the medium (usually 
air) at a velocity given by »v = B co. The velocity of light in vacuum is Co 
and the velocity of light in the medium is c = co/n where n is the index of 


refraction. The parameters used to describe the radiation are 


OF = Gicagien (2.3) 
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which is a constant, where yw is the permeability of the medium, and the 


radiation function R is given by 
iee— 277 sind I(u)P(k), (es) 


where @ is the angle between the direction of travel of the charged particle 
beam and the direetion of propagation of the emitted radiation, I(u) is ie 
diffraction function, and F(k) is the form factor. [Ref. 5: pp.2741,2742| 

The wave vector of the emitted radiation is k (k = w/c), and F(k) is the 
dimensionless form factor. That is, if p(r) is the charge distribution of a single 


bunch, then the Fourier components of the charge are 


p(k) = qF(k). (2.6) 


For a point charge, F(k) is identically one for all values of k. 

The dimensionless beam length parameter 7 which is defined as the ratio of 
the path length of the charged particle beam to the wavelength in the medium 
of the emitted radiation, appears in the radiation function explicitly, and 


implicitly as part of the diffraction parameter defined below: 


7 = L/ A. (ei 


1d 


The diffraction function is 


I(u) = sin u/u, (2.8) 


where the dimensionless diffraction parameter u depends upon the Cerenkov 
angle given by cos 0. = (n@)"', as well as upon the beam length parameter 


and the radiation emission angle: 


= 7) | (pe coseale (2.9) 


The parameter 7 determines the periodicity of the diffraction pattern and also 
determines its magnitude. Since R depends directly upon 7, the strength of the 
radiation is proportional to the square of the path length of the beam. [Ref. 5: 


pp. 2741,2742] 


D. FINITE-LENGTH PATH 

If the emission region has finite length, which may be realized by passing 
the electron beam through a gas cell, the radiation propagation direction is not 
confined to a sharp Cerenkov angle 4, but is spread over a range of emission 
angles. For the case of gas, for which the ordinary Cerenkov angle is a few 
degrees, the broadening is asymmetric about the Cerenkov cone, and the 
radiated power displays various interference lobes. This spreading effect depends 
only on the path length in the gas and does not depend on the electron bunch 
structure or periodicity of the bunches. [Ref. 3: pp. 3246,3251| 

Radiation from the polarized molecules along the path of the electron 
beam is formally equivalent to diffraction from a single slit at angles far from 


normal incidence. The shift in the radiation peak can be understood as a 
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diffraction effect arising from the linearly varying phase of the radiation along 
the beam interaction length [Ref. 3: p. 3246]. 

Considerable power (about half in the case of a gas) appears outside the 
main lobe, and the latter is the peaked at angles much larger than the 
Cerenkov angle. The Cerenkov power radiated per unit path length is increased 
by large factors approaching 2 orders of magnitude for a gas path of finite 
length compared to one of infinite length (Ref. 3: p. 3251]. 

This increase in power is associated with the finite length of the radiating 
medium only. The effect should be approximately the same for either single or 
periodic electron bunches and should occur whether the bunches are effectively 
point charges or have significant size. 

The value of the parameter u is affected by 7 and the beam energy 
through @. As either of these parameters is varied, the angular position of the 
zeros of I(u) changes. The emission threshold is defined as the situation when 
the first zero of the diffraction function enters the physical range. Applying 
this criterion, the result is that for small values of 7, a Cerenkov-—like 
radiation pattern is produced even when the product nf < 1. That is, the 
threshold occurs for n@ < 1 at small values of 7. As 7 increases (at constant 
nZ < 1), the threshold is crossed and the radiation pattern becomes more 


chaotic. A full discussion of these geometrical effects is given in Ref. 5. 


E. PERIODIC ELECTRON BUNCHES 

Cerenkov radiation was considered for periodic bunches of electrons such 
as would be emitted by a typical traveling-wave electron linear accelerator 
(LINAC). Because the distribution of intensity of Cerenkov is proportional to 


the frequency, the radiation at microwave frequencies would be low unless 


ly 


beams are intense and bunched so that coherent radiation by many electrons 
contributes (Ref. 2: p. 1531]. 

Periodic bunches of electrons produce Cerenkov radiation at harmonics of 
the bunch frequency, while a single bunch will radiate with a continuous 
frequency distribution. The effect of bunch size is reflected in the factor of the 
Fourier transform of the charge in the bunch, which causes the radiation to 
fall off at high frequencies, while at low frequencies all charges in the bunch 
radiate in phase. 

At low frequencies such that the wavelength of the e:nicted radiation is of 
the order of the bunch size, the electrons in the bunch radiate coherently. This 
leads to large enhancement factors for typical bunches consisting of 10° 
electrons, and in fact allows the radiation to be significant at microwave 
frequencies. [Ref. 3: p. 3246] 

At higher frequencies or harmonics, destructive interference, described by 
the Fourier transform of the charge density,decreases intensities with increasing 
frequency, until incoherent radiation takes over when the wavelength of the 
radiation is much less than electron spacing [Ref. 3: p. 3246]. 

When a single electron passes through a dielectric medium where a 
spatially periodic term is added to the dielectric constant, the result is 
radiation occurring even for electrons which do not exceed the velocity of 
light in the medium and at angles other than the Cerenkov cone angle. The 
non—Cerenkov part of the radiation is attributed to transition radiation. If the 
electron velocity were lower so that v/c were close to but less than unity, the 
peak in I would be pushed to the left toward 0° of @ such that cos 6 = v/c 


would be larger than 1. But the tails of the diffraction function I could extend 
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into the physical range 1 = cos 9 = —1, and this would be called transition 
radiation and be ascribed to the passage of the electrons through the 


boundaries of the gas cell. (Ref. 2: p. 1535] 


F. PROFILE DETERMINATION OF A ELECTRON BEAM BUNCH 
1. Effect of The Form Factor 

The bunch form factors are Fourier Transforms F(k) which differ 
depending on the charge distribution within a bunch. The form factor is 
identically one for a point charge, and for a finite charge distribution F(k) = 
1 for k = 0. If all the charge within a bunch is of the same sign F(k) must 
fall off as a function of k near the origin. Here we consider only line charge 
distribution and consequently the form factors are functions of the product of 
k, = kcos @ and appropriate lengths which are characteristic of the charge 
distribution of a bunch. 

Thus the behavior of the form factor F over the range of emission angles 
is dependant on the ratio of the charge bunch lengths to the wavelength of 
the emitted radiation. For a small value of these ratios, the value of F is 
essentially unity for all angles. Then the emitted radiation pattern is 
modulated only by the geometrical envelope function as described in Ref. 5. 

However for larger values of the ratios, the value of F over the range of 
emission angles can be significantly less than one. In this case the dependence 
of the form factor on angle is dramatic since the cosine term goes to zero at 
ninety degrees causing a form factor of unity. If, at same time the emission 
length ratio is small, a condition could be achieved where the majority of the 


Cerenkov power is radiated in a direction normal to the electron beam path. 
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A complete angular map of Cerenkov radiation patterns may be used to 
characterize the properties of relativistic beam pulses, and may enable one to 
determine the beam pulse charge distribution from measured radiated patterns 
[Ref. 4: p. 1996]. 

2 Bunch Shapes and Corresponding Form Factors 

Figure 2.3 diagrams the beam pulse shapes and presents the 


corresponding form factor expressions. [Ref. 8: pp. 16,17] 


Level function 





Sin kod 


kod 


Figure 2.3 Electron Beam Charge Bunch Shapes and 
Corresponding Form Factors. 
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Figure 2.3 (cont.) Electron Beam Charge Bunch Shapes and 
Corresponding Form Factors. 
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Figure 2.3 (cont.) 


Electron Beam Charge Bunch Shapes and 
Corresponding Form Factors. 
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where N = number of humps. 


Figure 2.3 (cont.) Electron Beam Charge Bunch Shapes and 
Corresponding Form Factors. 
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Ill. CALCULATION AND RESULTS 


A. PARAMETERS FOR CALCULATION 

Professor J. R. Neighbours and others of the Department of Physics, 
Naval Postgraduate School, Monterey, California, developed a computer 
program to calculate and plot angular maps of the radiated energy per unit 
solid angle for a given form factor. The program originally was designed to 
run on a Tektronix Computer and considered the following one dimensional 
(ie. in Z only) beam charge distributions: Gaussian, level (i.e. boxcar),level 
plus a sinusoidal, ripple combination,and a double hump. 

This work was aimed to modify the program to run on the IBM 370 
mainframe computer and plot the calculated power per solid angle, radiated at 
the frequency v, using the DISSPLA V10.5. The results are limited to the 
three form factors for the rectangular, Gaussian and trapezoidal charge 
distributions. Calculations for other charge distributions are left for future 
work. 

To calculate the total coherent power per unit solid angle, radiated at the 
frequency v by a periodic charged particle beam in traveling a finite distance 
L at constant velocity, a set of parameters were chosen close to those 
available in an experiment. 

(1) Index of refraction, n = 1.000268. 

(2) Speed of light in vacuum, co = 2.997925 x 10° m. 
(3) Fundamental frequency, v = 2.85 GHz. 

(4) Bunch charge, B = 1.5 x 107 Coulomb. 


These values give an n@ product of n@ = 0.98, which is 2 % below the 
infinite path length threshold value of unity. By choosing different radiation 
harmonics and varying the path length the radiation pattern above ( 7 = 12 ) 
and below ( 7 = 37 ) the threshold for 2.555 Mev electron bunches, could be 
investigated. The calculations were carried out at constant 7 for the various 


charge distribution as the size of the distribution was varied. 


B. GAUSSIAN FUNCTION 

Fig. 3.1 through Fig. 3.8 show the radiation intensity pattern of Gaussian 
line charge distributions in sub—Cerenkov radiation (n@ = 0.98 ) for two 
different value of 7. Fig. 3.1 through Fig. 3.4 are the results of using 7 = 12, 
with N = 3, varying the bunch length from 0.24 cm to 2.4 cm. The rest of 
them are for the case of using 7 = 37, and likewise, varying the bunch length. 
BL is the symbol used for bunchlength. 

Figures 3.1-3.4 show the effect on the emitted radiation pattern as the 
bunch length is increased while the beam path parameter 7 is held constant at 
n = 12. With a bunch length of 0.24 cm the radiation pattern is essentially 
that of a periodic set of point charges. As the bunch length increases, a hump 
first becomes evident in Fig. 3.2 when the wavelength of the emitted radiation 
begins to become comparable to the size of the bunch. Further increase in 
bunch size shows an increasing dominance of the form factor until finally the 
radiation pattern is almost completely dependant on the form factor. In Fig. 
3.4 the envelope of the radiation pattern is the form factor for the Gaussian 
charge distribution except for a slight skewing caused by the geometrical 


envelope function. 


29 


0.010 |; 
0.009 

0.008 | 
0.007 
0.006 : 


ONelels ee 


W (VVatt) 


0.C04 |} 
GO OGaa 
0.002 


GEO1S2) 





ee 


0 20 


igo! 


sub-Cerenkov 


E = 2.555 MeV 
7) = 2 

L = 42 cm 

A = 3.50 cm 
BL = 0.24 cm 


| 
Ht NWWWanne 


= ae -e 


40 60 80 100 120 140 160 
Angle (0, degrees) 


Gaussian Function with BL/A = 0.069, n = 12 


26 


180 


sub-Cerenkov 
E = 2.555 MeV 


7) 
L 
x 
BL 


9.0 | 


CmolO melas 
le con 
140 160 180 


LO CmN ZO 


Angle (@, degrees) 


| 
| 
| 
| 
BO 


. 
ee 
a 
i 
: 
f 
2 $ 
8 : 
¢ 
, ? 
tf ? t i 
‘i . 
¥ cai } 
* 
{ one ry ? 
, 
ae ry 7 
at a ‘ 
r | ‘ 1 
IS i 
i 4 ee tt 2 
p i ry 
ye d ; f 
LS a ee ae 
q 
$ Peer a We 
b ayy b 
y 4, | oi 
a ate 
) ry 
ae . bat S| 
j “{ f 
i a ro] 
i } 
| : If 
* : 


—— a ETDS CT TTL, CTT LS: A EE EE A, LT 


1 See Pee Se 
aa pe SS ET TY Ee TS ACT TT, SI 


Fig. 3.2 Gaussian Function with BL/A = 0.343, 7 = 12 





oO © O O O Oo Oo OQ Oo 
™ © LO + 4 rot = O 


(HIEAA) MA 


= 3. OOneina 
= 1.76 cm 


E = 2.599 VieV 
x 
BL 


1 = ae 


L = 4eu 


9.0 
oO 


7.0 





6.0 | 


120 140 160 180 


a ee Ly 
eee eee 


| 
| 
| 
i 


oS Bove 

EN The 

——— Try eee aa 
——> 

A a TF = 


ag Nip peo 
eee = 5 4 


Angle (@, degrees) 


ee 


| 
40 


WAN 


O 


2 
Yo 


4.0 
3.0 } 
Pat®, 
1.0 


0.503, n = 12 


Fig. 3.3 Gaussian Function with BL/A 


#10 ° 


W (Watt) 


10.0 | 


9.0 


8.0 


7.0 


6.0 | 
5.0 | 


4.0 | 


SA, 


220 


poets el 


40 


sub-Cerenkov 
E = 2.555 MeV 
ig) = VEZ 

L = 42 cm 

A = 3.50 cm 
BL = 2.40 cm 


ti 


60 80 100 120 
Angle {(@, degrees) 





mii 


3 i 


140 


| 


160 


Fig. 3.4 Gaussian Function with BL/A = 0.686, 7 = 12 


29 


180 


W (Watt) 


0.010 


sub-Cerenkov 


aes E = 2.555 MeV 
1) a 

0.008 L = 129%59ema 
A = 3.50 cm 
BL = 0.24 cm 


0.007 


0.006 


0.005 


0.004 


0.003 


0.002 


07008 





O 20 40 60 86 100 120 140 160 
Angle (©, degrees) 


Fig. 3.5 Gaussian Function with BL/A = 0.069, 7 = 37 


30 


180 


sub-Cerenkov 
= 2.555 MeV 


E 


37 
L = 129.5 cm 


(ee 


S50 Cm 
= 1.20c 


N= 


mM 


Be 

















(HIEAA) M 


Angle (©, degrees) 





fi 


0.343, 7 = 3 


Fig. 3.6 Gaussian Function with BL/A 


31 


129.5) cin 
eyo 0) (olag 
1 Seem 


E = Z.000eiey 


sub-Cerenkov 
ig) = 3) 


L 
N 
BE 


9.0 
8.0 
1.0 
6.0 


10.0 | 


‘ 
ot 
CQ) 
ai 
ord 
-s4 © 
= 
a7 
oe Cy 
mes, a 
men 
wwe 0s so ae 
a a -D= 
Co tes SPS) AES 
a : ~~ 
owe (the 
pe! fae ** em, - 
~Te- “eta = =: 4-8 - 
i . a Fa 
NT oo) mee 8 =e O 
= - - oa ~ 
ae eee FO 8 5 
Ce ee ee, Se, ee 0 dene oO ome. 
erage... §- Mee ee &® wo: =~ oe en AP _ 
wee ERP KR eg nl «a eee ee 
eens 6s a 8 -_ = «-42 . oN, Sa PT - = a 
ieee a a” 2 . i ie — 
ES EF ee 9 -7 me my gt ge - Peso ca Batters ee ae 
REE HOPS SED — tees Daas o ~~~ oem fe sty aa 
iene 8 me ee I 8 ae me AO Te emma 
eumsnt ea ewe OP OH Nw ene be ee TO’ ota" 9 a 2-eeag asta 
ST? on TP ST Ean * = gime ie - . = ew o ee eS pe) ON. Rete, oe 
re EE = Se ee ee CH 68s pee Ot ee meh EP Qe ke 
aE 268 & VO ~ ee BUN RR eme oe - ae Mae TEE Pom - = 2 eee ew 
ensayo eee EE Pe - EO OD. ee een O 
eee oer. PT Ue ter. eee ee ee ae Se He PRE OS, 2 mG tee 
ween 2 Oat rw we ee ee ee el hd o— tate oe Oe ee oe 
ee Oe er ek PAR de me SUHE Gem 2 me or ee one co 
Pe her ose Pate PO tee TF Sy LIL Tk ee eee OO or 
ED ee ae eR meee a A eS eee ee 
Ee ee I OB se Cute ~ Oe’ A Cee - ee ee eee 
SS A one oe BR ee Se me me 
ee eli ae -_-— = wae oe OPE? 4 ee UE 
nee i Lg ree. - eee wm REM). AOS F ee Oe we 
= = — a me =» a 
“Se -- Se teed _ = ow —_—_ - o- —ws ee eee - 
ope waweney = eS a. foe ee eet P| 
PRL OT ares © Ct at ae ms 
Ee Cie al BY 5 oa ee ———— 
: aa = ee > ee 2 
eT ee TA . = 
OEE ie es eee 
IE ee eT Te <a) 
: - en =. 
ew eA a — 
ae eee = OS a es = 
eR SE Tee 2 See Oe i 
~weawee= =o <a ae —T 
beer ee 
ss Mean esen a eee = 
eee a Pee ee = = g = 
we ag ot =e 
= - -eeeme- ee ie 
= — some oie 
TE Te ee tn O 
ES — a 
EARP eS ee an = 
—_e wz 
EE ng meme ees ~y 
erat 2 ** 
<< =). 
- mre, Es 
OE ene 
te 
i ee 


' 
| 
20 


5.0 
4.0 
3.0 
2.0 

1.0 


(IEMA) M 


7 


0.503, 7 = 3 


— 
— 


Angle (©, degrees) 
ae 


Gaussian Function with BL/A 


Pigs oat 


a 


50.0 


o Ol 


~ sub-Cerenkov 


E = 2.555 MeV 


1 = 3/7 


45.0 


L = 129.5 cm 
A = 3.50 cm 
BL = 2.40 cm 


40.0 | 


wom Bf eK ne Da 
2s eer we Ay OP 
eee me 8 eee ~~ 
: ee eet Ke ae 
Pky ek TS 
= wom - “*- We OE ews 
ROR ee = = 
ST Mi ener in - owen. Steet Yee 
eee a, ee ORE EET rt. Ste NTA ea) ¥ ln ameremnne So NON 8O — 
—neaand i SS ee TW ee gee “twa T MF. Tt _ } ee eee 
NE Sie i OE eT tee wee, A NOSE ee mare 
AE Fe PP. UP. OO mm me eS ee oe ee = ee ees Be ee 
emerson ewe. - etme + WETICE DF Ee gO ON en omete me ered IF A Ee ee ee 
ee ee hemline old alt ited ir I ee eet oa) Se = ante wr!’ ie we ee er Span ees tcommaan oe mae 
ee OE GAA i ee me BWSR - ee KC 0 en ee —= PS oP wee et te ee ee 8 9 Oe ee may : 
oe en ot a Oe a ee ew bed Mt Tas or hereneed I ee er ee og SE es oy 
= “+. | 8 woe sre Sm ye Ae yytts) at waa ca cee oo a i an at are me) Ou’ Cae Ese SE e s — 
ee 7. CUE AS Be! 6 i Ea oe Mer eer a6 o-— ese eee Dae Dee a mn ht gg OE es Oe ee - ee ee 
a Nee enna a stares ist tg h fOr ware “es samee- * act ~me EE a ae a 
ion os J me rere FW LSE CDCR BT Site ee RK SS ES OM CE ER. orem 2 Sp ee tm ee emegrretator com 
pote eeine Sita US AE SNE Oe ee OF oy ee pepewene pe wt US en ks Pe ee ee ed co 
TSS PUES DE HN ETC Be ee oem ree ‘= 3-stV eo: a. --— Se ee ee ae ow 
Tae re eh ENS OO TS Heres SE = Te erin es nape eR PE A I ree 
ES AO eee Ee othe AEP 8 tA 8e me - Sse? Sy CO BW nee eee ee ee 
= = 5 TO eee ete. ws. Zi a 
oe ee ee OE ee Bl ae PO Ye Soe EN ST ee hel OP a ee ee ee ng 
ee er ewe RE! CQ a Cee we nO NER OE 6a. me eee eee 
TO en es tree PI ae Se Pe AE. he OR meres eres Os 
= . — —— = 
AE OME aera ste eens = 2S ae coe oe oo anes OS 
= we kg ale ne Be a AD 0 EL Nes arene 
pepe ee TAA Ere SD em er ere i Soe eaatepepeenl 
Se ee one EE Se amon a 
=~ =Sr RA eee Se +s TE. > 
e-Pea pe oo —- 22 Re. ewe ey = 
: => Fe ee RE ant lee 
een > A 
5 4 me ee. ee 
Pq oe 


So. 0g 


30.0 


25.0 |- 


(HIEAA) MA 


20.0 


S18) 


10.0 





0 


Angle (@, degrees) 
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Similar remarks can be made for the longer finite path results shown in 
Fig. 3.5 — 3.8 for 7 = 37. In this set of figures, because this value of 7 causes 
the radiation to be below threshold, the radiation patterns are wildly 
oscillating — almost chaotic. However exactly the same thing happens to the 
envelope of the radiation. Fig. 3.8, calculated for a bunch length of 2.4 cm, as 
is Fig. 3.4, also shows the radiation envelope to be the form factor for the 


Gaussian charge distribution. 


C. LEVEL FUNCTION 

Fig. 3.9 through Fig. 3.22 show the radiation intensity lobe development 
as the level function charge distribution length varies for the same two values 
of n. Fig. 3.9 through Fig. 3.15 are calculated for 7 = 12, the remainder 
through Fig. 3.22 for 7 = 37. LFL is the symbol for the level function length. 

In Fig. 3.9, where the level function length is much smaller than the 
radiation length, the radiation pattern is almost the same shape as Fig. 3.1 
which is for the Gaussian charge distribution. Proceeding through to Fig. 3.13, 
as the length of the charge distribution is increased the form factor becomes 
prominent. Fig. 3.13 is of interest in that it also shows considerable radiation 
at small angles. This result occurs because the form factor for the level 
function has secondary maximum which at low angles combine with the 
relatively large value of the diffraction function to give substantial radiation. 

Further increases in level function length give an increased role to the 
form factor until finally in Fig. 3.15, it has become dominant. Here the 
envelope of the radiation is seen to be the (skewed) form factor of the level 
function. The small lone peak at approximately ten degrees occurs at the 


maximum of the geometrical envelope function. 
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Fig. 3.20 Level Function with LFL/A 
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Fig. 3.16 to Fig. 3.22 drawn for 7 = 37 show the increasing dominance of 
the form factor below threshold. The pattern oscillates very rapidly, but 
otherwise the results are as above. In fact the rapid oscillations make the form 


factor more evident in this series of figures. 


D. TRAPEZOIDAL FUNCTION 

Both the Gaussian and level function distributions have a_ single 
characteristic length; the trapezoid is more complicated in having two, the 
length of the base, lp, and the length of the top, ly. In the fourier transform 
these appear as f = ( lp + It )/2 and 1 = ( lp -— 1k )/2. 

Fig. 3.23 through 3.32 show radiation patterns for a trapezoidal charge 
distribution for different values of top and base lengths. All except the last 
three, Fig. 3.30 through Fig. 3.32 are for 7 = 12. TFTL is the symbol for the 
top level length of the trapezoidal function and TFBL, for the base level 
length of trapezoidal function. 

As mentioned earlier, the form factor is sensitive to the ratios of the 
characteristic lengths to the wavelength of the emitted radiation. Here the 
appropriate ratios are ¢/n and 1/A. In Fig. 3.23 both these ratios are less than 
one so that the radiation pattern is characteristic of a point charge 
distribution. Fig. 3.24 drawn for ¢@/A 1.2 and 1/A = 0.29 shows a different 
pattern with the form factor already becoming quite evident. 

Fig. 3.25 through Fig. 3.27, calculated for a constant slightly larger ratio 
of ¢/A = 1.71, and increasing values of 1/X are all similar since the radiation 
wavelength is too large to be sensitive to the smaller details of the charge 
distribution. The same comment applies also to Fig. 3.28 for which @/A is still 


larger but 1/X is less than one. 
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The last figure of the 7 = 12 series is drawn for both @/X and 1/A greater 
than one. Consequently one would expect the radiation envelope to display 
both terms in the form factors for a trapezoidal charge distribution. That it 
apparently does not may be due to the fact that the lengths chosen in order 
to achieve these make for a charge distribution which is almost triangular. 

Figures 3.30-3.32 are calculated for the below threshold condition of 7 = 
37. All have @/X greater than one and 1/A less than one. As with other figures 
with this value of 7, the envelope is more evident. 

When both values of @/X and 1/A increase the peak values decrease mainly 
in the region of angles between 0° and approximately 55° and the radiation 
intensity is enhanced around 90°. It is difficult to determine the development 
of the radiation pattern with respect to ¢/A and 1/A because both values can 
change the radiation pattern at same time. To analyze the radiation pattern in 
another way the angle where zeros happen can be investigated by calculation. 
Let's set: 


Kel ong, Se (3.1) 


where kz = k cos @, and m,n are integers. Substituting for k gives 
cos = #n—, cos fy = kM (3.25 


If aS in most cases considered here < 1, these are no solutions for @p. 
xe 
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Thus only the solutions for #, are considered below. A example is Fig. 3.25 


where <* = | ee —sOo Cine — o-o CM. 





S5) 1 
cos @ = £#n Tae = tn Tawa = +n x 0.58 foe3)) 


The solutions are for @, are 
f - 545° and 64 = 125.5°. (3.4) 


which agree with the figure. These values also apply to Fig. 3.26, Fig. 3.27, 
Fig. 3.31, and Fig. 3.32 which have the same ratio of 
In Fig. 3.29 where —~ 1.29, I/A = 1.27, 
ot — a cos @ m= (20) 
. 20 e2 an Wee 
Miessoluitions are 6; = 39.2°, 6, = 129.2° for @, and @, = 38°, @., = 128°. In 


This case the zeros for 6, and 9, are so close together as to be 


indistinguishable in the figure. 
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IV. CONCLUSIONS AND RECOMMENDATIONS 


The object of this work is to assess the effects of the charge distribution 
of a bunch on the sub-Cerenkov radiation from a charged particle beam, 
specifically investigating the short beam path region at energies above and 
below the threshold. It contains many calculated radiation patterns for the 
three kinds of charge distributions. 

With the charge distribution of Gaussian function, the radiation intensity 
envelope near 90° represents the square of the bunch form factor F(k) when 
the bunch length is larger than a half of the wavelength. When the bunch 
length is about same as the wavelength or larger than that of the level charge 
distribution the radiation pattern around 90° is also dominated by the square 
of the bunch form factor F(k). 

It is difficult to determine the characteristics of the radiation pattern on 
the trapezoidal function by the ratio of bunch length to wavelength because 
the level top length and level base length can change the shape of radiation 
pattern at same time. However the radiation pattern of the trapezoidal 
function can be partially characterized by analyzing the zeros of the curves. 

A summary of results is : 

le When the bunch length is much smaller than the wavelength of the 
emitted radiation, the radiation pattern is characteristic of point 
charge. 

2. When 7 is increased the total radiation pattern does not change 


except for minor differences in the frequency of oscillation. 


ay When the wavelength of the emitted radiation is comparable to or 
larger than the bunch length, a significant amount of radiation 
appears at near 90° to the beam. 

4, The envelope of the 90° radiation is the square of the bunch form 
factor. This envelope is skewed slightly as a result of monotonic character of 
the geometrical envelope function. | 

To enhance this preliminary work the following is recommanded: 

1. Investigate the behavior of the square of form factor F(k) only when 
the bunch length varies. This will give the insight into the behavior 
of radiation pattern. 

2. Determine why varying the top length of the trapezoidal pulse shape 
causes significant enhancement of the forward lobe envelope. 

3: Determine the effects of changes in beam energy, pulse width,pulse 
frequency, and path length. 

4. Confirm the results found by mapping the microwave radiation 
measured at a comparable particle accelerator. 

The examples presented here are idealized and there has been no attempt to 
take into consideration such factors as beam pulse erosion as it transverses the 
atmosphere and the ground plane reflection in interference[Ref. 4: p.1996]. But 
the complete angular map of the envelope of the radiation pattern could 
provide a measure of F(kz) and an excellent reference to establish a method of 
determining the beam pulse charge distribution based on sub—Cerenkov 


radiation patterns. 


63 


APPENDIX A 


PROGRAM FOR CALCULATION OF RADIATED POWER 


He Ht 4 4 HE HE 3G HE HE HEE HE EE HE HE HE 4G HE 9G 9 FE HE HE HE HE HE 3 EE HE HE EE FE HE HE EE HE HE EE HE HE FE FE HE HE HE HE HE HE HE HE EE HE HE EE FE HE OE OE EE EE He 


* THIS PROGRAM IS TO CALCULATE THE POWER PER UNIT SOLID ANGLE RADIA- * 


* TED AT THE FREQUENCY NU. IT WAS ORIGINALLY DEVELOPED BY PROFESSOR 
* JOHN R. NEIGHBOURS AND OTHERS FOR TECTRONIX COMPUTER. THIS AUTHOR 
* MODIFIED TO RUN ON IBM 370 MAINFRAME COMPUTER IN DEC. 1987. 

* THE THEORY AND ALL EQUATION ARE WELL GIVEN IN REF. 5. 


* * *k * 


Fe ME HE HE EME HE HE HE HE HE HE HE EE 9b EE HE HE HE EE HE HE IE FE HE HE 4 HE EE HE EE 3b FE HE EE 3 HE HE HE HE HE FE FE HE HE HE EE HE HE HE EE HE EE HE HE 0 EE EE 


PROGRAM JNN 
/VA/HH,HU,SDXX /DOUPA/L,ETA,TMAX,DT,E,N,B,LL,RR,NN,RFTL, 


COMMON 

i: 
DOUBLE 
DOUBLE 
DOUBLE 
DOUBLE 
DOUBLE 
DOUBLE 
DOUBLE 
DOUBLE 
DOUBLE 
REAL 

+ 


RFBL,TFIL, TFBL, TRIBE Nees 


PRECISION C,CO,MU®,NU0,0,P1,50,8,NBETA, 

PRECISION TMAX,DI,NU,L, ETA ,N, be, Like NN 

PRECISION IIMAX,KK,CONST,WMAX,THETA,W,T,KKZ,H,DENOM1, 
PRECISION DENOM2,K®@ 

PRECISION GAMMA,V,U,1I1,R,HA,HB,HC,HD,DD,DEXxIN 
PRECISION REFTL,RYSL, TPL, Deelah eee 

PRECISION DENOM3,DENOM4 , DENOMS , DENOM6 , DENOM7,DENOM8, 
PRECISION Soh he HGH t 

PRECISION HJ,HK,HL,Z,LI,LIN,HA AW ke Sp 


AWMAX , XTHETA, XW, YMAX,SDYY,SDY ,JX,3YY (YY , YUM Xe 
XX ,XNUM,RTMAX,RSDXX,RWMAX,RTHETA, RW 


REAL REL,REE,REN,REB,RELL,RERR,RENN, DELTA 
REAL RERFTL,RERFBL, RETPTL, RETPSBL, RETRIS, ReaD 


Q 


Q 


INTEGER 
CHARACTER 


I,K,KMAX, IMAX,HH,J,JXMAX, JYMAX,IPA,P,HU 
*1 PARACH,EXPPTS, EXPTCH,SAVEQ,PPRE,CHANGE,ANS,AXCH, 
ANS2 , ANS3 


CHARACTER*13 NAMEPF ,NAMEWF , NAMEEX 
DIMENSION XTHETA(1:100) 
DIMENSION XW(1:1900) 

DIMENSION THETA(@:45900) 
DIMENSION W(@:4500) 

DIMENSION H(1:7) 


* ASSIGNMENT OF CONSTANTS 


C=2.997 


122D08 


CO=2 .997925D08 


MU@=1.2 


PIl=45.14 


566437061D-96 


HH HHH HH MH HH TO VARY CHARGE HH HE EH HE HE HE He HE 


159265556979 5200 
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BO=0.5117D@ 


2D EN ae ene 
PRINT*, ' ENTER BUNCH CHARGE.(DEFAULT= 1.5E-12) ' 
READ*, Q 
IF(Q.GT.@) GO TO 45 
Q=1.5D-12 
45 PRINT*, ' ENTER FUNDAMENTAL FREQUNCY NU@. (DEFAULT= 2.8559)’ 
READ*, NU@ 
PECNUGO.GT.0) GO TO 55 
NU@=2 .85D89 
2D CONTINUE 


PRINT*,’THE INDEX OF REFRACTION IN THE AIR IS 1.000268(DEFAULT).’@ 
PRINT*,’DO YOU WISH TO CHANGE? (Y/N) ' 
READ ‘(1A)’, ANS2 
IF(ANS2.EQ.’'Y’) THEN 
PRINT*,’ENTER THEINDEX OF REFRACTION. '° 
READ*, DEXIN 
ELSE 
DEXIN=1.808268D0 
END IF 
C=C@/DEXIN 


C ieee DEGESTON OP PARAMETER FILE*****«** 


PRON, | * 

PRN *, 1: CREATE NEW PARAMETER FILB’ 
PEN 2 eUSHsOLD FLLs. 

PRON",  ' O28 salty GCIs: 

PRONT * 

PEN, ENTER YOUR CHOICE <:’ 

READ '(1A)’, PARACH 

IF(PARACH.EQ.’3') GOTO 700 

IF(PARACH.EQ.’2’) CALL INAMPF 
IF(PARACH.EQ.'1’') CALL MAKEPF 


* 


C 


Se 


- READ IN ALL PARAMETERS FROM EXTERNAL me 
- FILE/PRINTOUT PARAMETERS ia 


HH Ht Ht FE Et OE OE OE HE OE 3 OE EE Ee ee HH HH HHH 

c 

99 WRITE( 6,100) 

100 FORMAT( ///’THE PARAMETERS YOU ENTERED ARE:’ ) 
REWIND 1 

120 READ(1,*) HH 

Sl WRITE(6,135) HH 

155 FORMAT(/’FOURIER TRANSFORM CHOICE : ieee) 
IF (HH.NE.1) GOTO 137 
WRITE (6,136) 

136 FORMAT( /,37X,’GAUSSIAN FUNCTION ' ) 

137 IF (HH.NE.2) GOTO 139 
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12S 
159 


WRITE (6,138) 


FORMAT (/,37X,'LEVEL FUNCTION 


IF (HH. NE <5) GOLOI4) 
WRITE (6,140) 


FORMAT (/'SUM OF LEVEL AND RIPPLE FUNCTIONS’ ) 


IF (HH.NE.4) GOTO 143 
WRITE (6,142) 


FORMAT(/’MULTIPLE HUMP FUNCTION’ 


IF (HH. NE. 5) 3GOLe mia 
WRITE (6,144) 

FORMAT (/'ROUNDED FUNCTION' 
IF (HH. NE 36) GOTOmI=7 
WRITE (6,146) 


FORMAT (/'TRAPEZOIDAL FUNCTION’ 


IF (HH.NE.7) GOTO 149 
WRITE (6,148) 
FORMAT (/’TRIANGLE FUNCTION’ 


READ(1,*) ETA 

WRITE(6,151) ETA 

FORMAT( /' ETA (L/RAMMDA ) 
READ(1,*) TMAX 
WRITE(6,155) TMAX 


FORMAT( /' MAXIMUM THETA VALUE 


READ(1,*) DT 
WRITE( 6,165) DT 


FORMAT( /’ THETA CALCULATION INCREMENT 


READ(1,*) SDXX 
WRITE(6,175) SDXX 


FORMAT( /’ THETA MARKING INCREMENT 


READ(1,*) E 

READ(1,*) NBETA 
WRITE(6,185) E 
WRITE(6,185) NBETA 

FORMAT( /? BEAM ENERGY 
FORMAT(/’ VALUE OF NBETA 
READ(1,*) N 
WRITE(6,195) N 

FORMAT( /? HARMONIC NUMBER 


IF CoE 1) GOLO=206 
READ(1,*) B 
WRITE(6,2905) B 


FORMAT( /? BUNCH LENGTH (M) 


IF CHE bt 62) GO LOg 216 
IF (HE -GT 4) SGOreeziG 


READ(1,*) LL 
WRETEC6,215) 05 


) 


) 
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2 
e 


, Die) 


»D16.4) 


»D16.4) 


(16.4 > 


' D1644) 


Dito 


5 D16.43) 


',D16.4) 


FORMAT ( /' LEVEL FUNCTION LENGTH (M) : ',D1i6.4) 


i CH ONE. 3) COTO19255 


READ(1,*) RR 
WRITE(6,225) RR 
FORMAT(/’ FOURIER RATIO CONSTANT : D16.4) 


READ(1,*) NN 
WRITE(6,235) NN 
FORMAT(/’ NUMBER OF RIPPLES : P; DiG4s) 


IF (HH.NE.4) GOTO 236 


READ (1,*)HU 
WRITE (6,9239) HU 
FORMAT( /’ NUMBER OF HUMPS: ie 


Peete Ep oy eee 24 | 


READE C1," ) REIL 
WRITE (6,238) RFTL 
FORMAT (/’ ROUNDED FUNCTION LEVEL TOP LENGTH 


READ (1,*) RFBL 
WRITE (6,240) RFBL 
FORMAT (/’ ROUNDED FUNCTION BASE LENGTH 


IF (HH.NE.6) GOTO 246 


READ C1, * ) TFTL 
WRITE (6,243 )TFTL 
FORMAT (/' TRAPEZOIDAL FUNCTION TOP LEVEL LENGTH 


READ (15%) TEBL 
WRITE (6,245)TFBL 
FORMAT (/' TRAPEZOIDAL FUNCTION BASE LENGTH 


Tesh ONE. ee GOTO 249 


READ (1,*) TRIBL 

WRITE (6,248)TRIBL 

FORMAT (/? TRIANGLE FUNCTION BASE LENGTH 
WRITE(6,25@) 


FORMAT(//'DO YOU WISH TO CHANGE ANY OF THESE ?(Y OR 


READ’(A1)’, CHANGE 


IF(CHANGE.EQ.'Y’) GO TO 700 


HH HH HE HE HH IE EH HE HE EH He HH HH HH HH HH 


eH HEE 


BEGIN CALCULATION *** * 
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DON 


2eeDi1.6 


Do 


Dio 


= Dien 


4) 


4) 


a) 


wa) 


4) 


LM i a a ee ee D 
C 
p22 ITIMAX=TMAX/DT 

IMAX=INT( IIMAX ) 


NU=N*NUQ 

IF ((HH.GE.2).AND.(HH.LE.4)) THEN 
DD=LL/2.@D0 

END IF 

é 

KK=2.0DQ*PI*NU/C 

CONST=MUG*C* (NUQ**2.0)*(Q**2.0)/(8.0DO*(PI**2.0)) 

WMAX=0.90D@ 

DO 350 I=0,IMAX 
THETA(I)=DBLE(I)*DT 
T=PI*THETA(1I)/180.0D0 
KKZ=KK*DCOS(T) 

é 
IF (HH.NE.1) GOTO 326 
C 
H(1)=DEXP(-1.@D0*((KKZ*B/2.0D0)**2.0)) 
€ 
326 CONTINUE 
IF (HH.LT.2) GOTO 328 
IF (HH.GT.3) GOTO 328 
C 
DENOM1 =KKZ*DD 
IF (DENOM1.EQ.0) GOTO 327 
H(2)=(DSIN (KKZ*DD) )/(KKZ*DD) 
GOTO 328 
527 H(2)=0.90D0 
Cc 
328 IF(HH.NE.3) GOTO 333 
C 
KO=PI*NN/DD 
DENOM2=(KKZ-KO )*DD 
€ 
IF(DENOM2.EQ.@) GOTO 329 
C 
HA= DSIN(DENOM2 ) /DENOM2 
C 
GOTO 330 
c 
329 HA=0.90D@ 
330 DENOM3=(KKZ+KO)*DD 
C 
IF (DENOM3.EQ.0) GOTO 331 
Cc 
HB=DSIN(DENOM3 )/DENOM3 
6 
GOTO 332 
C 
331 HB=0.0D0 
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Bo2 H(3)=H(2)+((RR/2.@D@)*( HA+HB) ) 


C 
po ipa CHH NE 4 eGOL@ 555 
C 

Z=(KKZ*DD) 
C 

IF(Z.EQ.@) GOTO 334 
C 


H(4)=1.0D0*((DSIN(Z))/Z) 

DENOM4= 1.0 

DO 9333 P=1,HU 
KE=(2.@D@*DBLE(P))-1.0D0 
LI=((2.@D@*DBLE( HU) )-KE)/((2.@D@*DBLE(HU))-1.@D0) 
LIN=((-1.0D®)**DBLE(P ))*(-1.0D@) 
HX=(LIN*LI*DSIN(Z*LI))/(Z*LI) 
HW=LIN*LI 
H(4)=H(4)+HX 
DENOM4= DENOM4+HW 

5555 CONTINUE 
H(4)=H(4)/DENOM4 


C 
GOTO 335 

c 

334 H(4)=0.0D0 

@ 

555 IF (HH.NE.5) GO TO 339 

c 
HL=(8D0*16D@)/(((RFBL-RFTL)**2.0)*(KKZ**3)*(RFBL+RFTL) ) 
HG=HL*DSIN( (KKZ/4.9D@)*(RFBL+RFTL ) ) 
HJ= (RFTL*DSIN( KKZ*RFTL/2.0D@))/(KKZ*RFTL/2.@D0) 
HK= (RFBL*DSIN(KKZ*RFBL/2.0D@))/(KKZ*RFBL/2.0D0) 
HI=(32D0/((((RFBL-RFTL)*(KKZ))**2.0)*(RFTL+RFBL) ))*(HJ+HK) 
H(5 )=HG-HI 

C 

339 IF (HH.NE.6) GOTO 345 

C 
DENOM6=KKZ*(TFBL+TFTL )/4.@D0 
DENOM7=KKZ*(TFBL-TFTL)/4.0D@ 
HE=DSIN( DENOM6 ) / DENOM6 
HF =DSIN(DENOM7 ) /DENOM7 
H(6)=(TFBL+TFTL )/2.0D0*HE*HF 

C 

345 IF(HH.NE.7) GOTO 349 

@ 
DENOM5 =(KKZ*TRIBL)/4.90D@ 

C 
IF (DENOM5.EQ.@) GOTO 346 
H(7) = (DSIN(DENOMS )/(DENOMS ))**2.@ 

C 
GOTO 349 

C 

346 H(7)=0.0D@ 
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349 CONTINUE 
C 
C GAMMA =E/E@ 
L=ETA*C/NU 
C GAMMA=DSQRT(1.@D0/(1.0D0-(NBETA/DEXIN)**2) ) 
C V=C@*( DSQRT(1.9DO-(1.9DG/(GAMMA**2.0)))) 


V=CO*(NBETA/DEXIN) 
U=(KK*L/2.0D0)*((C/V)-DCOS(T) ) 
II=DSIN(U)/U 
R=KK*L*H(HH)*II*DSIN(T) 

W( I )=CONST*(R**2.0) 


C 
C *** CALCULATION OF Serum y= = 
C 


Ate WMAX=DMAX1(WMAX, W(T) ) 
C 

350 CONTINUE 

C 


C eee DECISION TO SAVE VALUES CALCULATED ****i0= 


WRITE( 6,360) 
360 FORMAT(/'DO WANT TO SAVE THE RESULTS IN A FILE? (Y/N)’ ) 
READ'( Ai)’, SAVEQ 


IF(SAVEQ.EQ.'N') GOTO 440 


PRINT *, ‘UNDER WHAT FILE NAME? ' 
READ '(A10)’, NAMEWF 
OPEN( UNIT=3,FILE=NAMEWF,STATUS='NEW’ , FORM='FORMATTED’ ) 
DO 370 I=0,IMAX 
RTHETA=REAL(THETA(I)) 
RW=REAL(W(I)) 
WRITE(3,380) RTHETA,RW 
380 FORMAT(4X,F8.3,2X,F20.18) 
370 CONTINUE 
WRITE( 3,385) 
385  FORMAT(7X,'@.000' ) 
WRITE(3,390) WMAX 


39@ FORMAT(2X,'WMAX = ',Di6.4) 
WRITE CS, 292 Jeera 

392 HORMATC2 KE UA = 2". DiG.4)) 
WRITE(3,394) NBETA 

394 FORMAT(2X,'NBETA = ’,D16.4) 

C WRITE(3,396) B 

C396 FORMAT(2X,' BUNCH LENGTH = ',Di6.4) 

C 


C weeeneee = DECISION TO PRINTOUT,RE-RUN,OR EXD? 5 
C 


440 WRITE(6,445) 
445 FORMAT(////? 1: PRINT OUT VALUES’ ,/ 
+ 2: RUN PROGRAM AGAIN’, / 


70 


+ 3: EXIT’,/// 
+ ’ ENTER CHOICE. ') 
READ’(A1)’, PPRE 


C 
PRCPPRE- BO 1) GOloO 660 
PRCPERE SEG. 12°) GOTO 05 
PECEPRE BO. °S’ ) “GOTO 702 
GOTO 440 
C 
fee 2 PRINT OUT VALUES *** 
C 
660 WRITE( 6,670) 
670 FORMAT(’ THETA WeefeaQ7=")/) 
DO 690 I=0,IMAX 
WRITE(6,680) THETA(I),W(I) 
680 FORMAT(F8.2,F15.4) 
699 CONTINUE 
C 
'o HHRHHKRHH HH EH PROGRAM BNDS HHH HH HH HH HH HH HH HH HH HH HH HH HHH HHH 
700 END 
C 
C 
C 


Ht 4G HE HE HE FE HE FE TE HE HE HE HE HE HE HE HE HE EH HH 


* *** SUBROUTINE INAMPF *** * 
He HH HH HH HE HE HE OE He Ee OE OE HE OH HE HO HH He HH 
C 
SUBROUTINE INAMPF 
c 
C *** INPUT OLD PARAMETER FILE NAME 
CHARACTER*10 NAMEPF 
WRITE( 6,990) 
900  FORMAT(/’ENTER THE NAME OF THE PARAMETER FILE YOU WISH TO’, 
+ ’ USE(MAX 10 CHAR) ° ) 
READ’(A10)’, NAMEPF 
WRITE(6,911) 
911 FORMAT(/’TEST LINE ONE’ ) 
OPEN( UNIT=1,FILE=NAMEPF,STATUS=’OLD’ ,FORM='FORMATTED’ , 
7 BLANK='NULL’ ) 
WRITE(6,912) 
912 FORMAT(/’TEST LINE TWO’ ) 
RETURN 
END 
C 


At i 


fee 6 SUBROUTINE INAMEX *** * 


He HH HH HE HE HE HE OE HE HE HE HE EE HE HE HHH 


SUBROUTINE INAMEX 
C 


C *** INPUT OLD EXPERIMENTAL POINTS FILE NAME *** 
CHARACTER*1® NAMEEX 


al 


WRITE( 6,190) 
19820 FORMAT(/’ENTER NAME OF EXPERIMENTAL POINTS FILE YOU WISH’, 
+ ‘TO USE(MAX 18 CHAR) ' ) 
READ’(A1@)’, NAMEEX 
OPEN(UNIT=2, FILE=NAMEEX,STATUS=’OLD’ ,FORM=’FORMATTED’ , 
+ BLANK=’NULL’ ) 
RETURN 
END 
C 
C 


tHE HE HE HE HE HE HE EE HE HE EE OE HE EE OE HE HE OE OE OE OE Ot OE 


* ### SUBROUTINE MAKEPR *** * 


He HE HE HE HE HE HE HE HE HE HE HE HE HE HE HE HE HE OE HE HE OE Oe HE HE OH HE OH 


SUBROUTINE MAKEPF 
C 


C *** INPUT OF PARAMETERS INTO NEW FILE (UNIT=1)*** 


COMMON /VA/HH,HU,SDXX /DOUPA/L,ETA,TMAX,DT,E,N,B,LL,RR,NN,RFTL, 


+ RFBL,TFTL,TFBL,TRIBL,NBETA 
DOUBLE PRECISION L, ETA, IMAA,DI,&,N, 3, LE SDA ican. 
+ RR NN, RETO REESE, TP, IP Bia tater 


INTEGER HH, HU 
CHARACTER*1@ NAMEPF 
WRITE(6,13090) 
13909 FORMAT(/’ENTER THE NAME FOR THE NEW PARAMETER FILE.’, 
+ ‘(MAX 10 CHAR) ° ) 
READ’(A19)’, NAMEPF 
OPEN(UNIT=1,FILE=NAMEPF,STATUS=’NEW’ ,FORM=’FORMATTED'’, 
+ ACCESS=’ SEQUENTIAL’ , BLANK=’NULL’ ) 
F ieee : 
WRITE(6,13290) 
1320 FORMAT(/////’NOTE: ALL THE PARAMETERS MUST BE ENTERED’, 
+ *’ ALTHOUGH SOME MAY NOT BE USED’ ) 
WRITE(6,13390) 


1339 FORMAT(/’ 1: GAUSSIAN FUNCTION’/’ 2: LEVEL FUNCTION’, 
+ ‘ie 3: SUM OF LEVEL AND RIPPLE FUNCTIONS’, 
+ a 4: MULTIPLE HUMP FUNCTION’, 
+ (2 >: ROUNDED FUNCTION’, 
+ /’ 6: TRAPEZOID FUNCTION’, 
+ fe 7: TRIANGLE FUNCTION’ ) 


WRITE(6,134@) 
1349 FORMAT(/’ENTER THE FOURIERTRANSFORM CHOICE: ’ ) 


READ*, HH 

c 

1356 IF(HH.LE.7) GO TO 1359 

Cc 
WRITE(6,1357) 

1357 FORMAT(’ FOURIER TRANSFORM CHOICE NOT DEFINED.’///, 
- , ENTER NEW CHOICE:’ ) 
READ*, HH 

S 
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SO-10 1556 


'559 “PRITE(1,*) HH 
Wel TEC6, 1579) 
1370 FORMAT(/’ENTER ETA (L/RAMMDA). ( DEFAULT = 12.0D-2 ) ’ ) 
READ*, ETA 
Pelan Gi. O10) GOLO 1385 
ETA=12.@0D@ 
1385 WRITE(1,*) ETA 
WRITE(6,1400) 
1400 FORMAT(/’ENTER THE MAXIMUM VALUE OF THETA TO BE’, 
+ ' CALCULATED(DEG).'/’ ( DEFAULT = 45.0 ') 
READ*, TMAX 
IF(TMAX.GT.@.@) GOTO 1415 
TMAX=45 . ODO 
1415 WRITE(1,*) TMAX 
WRITE(6,143@) 
143@ FORMAT(/’ENTER THE ANGLE CALCULATION INCREMENT( DEG)’, 
+ TOE AUET=0.1) * ) 
REA, DT 
Pe DT.Cl.0.0) GOTO 1445 
DT=0.10D0 
eas WRETEC1,*) DT 
WRITE(6,1460) 
1460 FORMAT(/’ENTER THE ANGLE MARKING INCREMENT(DEG)’, 
+ OPER AULT=15.50) —’ ) 
READ*, SDXX 
LE CoDKX. CE.0.0) GOTO 1475 
SDXX=15.0D0 
1475 WRITE(1,*) SDXX 
WRITE(6,1490) 
C1490 FORMAT(/’ENTER THE BEAM ENERGY (MEV) (DEFAULT=100.0) ’ ) 
1490 FORMAT(/’ENTER THE VALUE OF NBETA (DEFAULT=@0.98) ’ ) 
C READ*, E 
READ*, NBETA 


C mer .CT.0:0) GOTO 1595 
IF(NBETA.GT.0.8) GOTO 1505 
C E=100.0D@ 
E=@.98D0 


e505 WRITE(1,*) E 
1505 WRITE(1,*) NBETA 
WRITE(6,15290) 
1520 FORMAT(/’ENTER THE HARMONIC NUMBER (DEFAULT=1.0) ’ ) 


READ*, N 
beCNEGT.0.0) GOTO 1535 
N=1.0D@ 
feoe WRITE(1,*) N 
C 
Poe enh | GOTO. 1575 
C 


WRITE(6,1550) 
1559 FORMAT(/’ENTER THE BUNCH LENGTH (M).( DEFAULT = @.24D-2 ) ' ) 


1) 


1680 


+ 


REAL; B 
IF(B.GT.9.90) GOTO 1565 


B=0.24D0 
WRITE(1,*) 8 


IF ((CHH.LY.2).CR (HH Cl 4) col Cmirewe 
WRITE(6,1580) 
FORMAT(/’ENTER THE LEVEL FUNCTION LENGTH (M).’, 
* (DERAULT=93019 5.) 
READ*, LL 
IF( LL. GT 080) GOTO 1595 


LL=0.01D0 
WRITE(1,*) LL 


i? (HE: NE. 5) "GOre@iesS 
WRITE(6,1610) 
FORMAT(////’ENTER THE FOURIER RATIO CONSTANT (DEFAULT=1.90).’) 
READ" 7 RA 


IFPCRR.GT 930 ) "GOO io25 


RR=1.0D0 
WRITE(1,*) RR 


IF (HH.NE.3) GOTO 1665 

WRITE(6,1640) 

FORMAT(/*ENTER THE NUMBER OF RIPPLES (DEFAULT=1.0) ° ) 
READ*, NN 


IF(NN.GT.@.9) GOTO 1655 


NN=1.0000 
WRITE(1,*) NN 


IF (HH.NE.4) GOTO 1700 

WRITE(6,16790) 

FORMAT(/’ENTER NUMBER OF HUMPS (DEFAULT=2) ° ) 
READ*, HU 


IF(HU.GT.2) GOTO 1680 


HU=2 
WRITE(1,*)HU 


74, 


FORMAT(/’ENTER THE ROUNDED FUNCTION LEVEL TOP LENGTH (CM)’, 


FORMAT (/’ENTER THE ROUNDED FUNCTION BASE LENGTH (CM).’'/, 


FORMAT (/’ENTER THE TRAPEZOIDAL FUNCTION LEVEL TOP’, 
(DEFAULT 


FORMAT (/’ENTER THE TRAPEZOIDAL FUNCTION BASE LENGTH (CM)’, 


FORMAT (/’ENTER TRAIANGLE FUNCTION BASE LENGTH (CM).'//, 


1.0) 


2.0) 


eyo) 


iD) 


IF (HH.NE.5) GOTO 1800 
WRITE (6,171@) 
+ /? (DEFAULT = 
READ*, RFTL 

IF (RFTL.GT.@.@) GOTO 1730 
RFTL = 1.0D0 

WRITE (1,*) RFTL 

WRITE (6,1760) 
n (DEFAULT = 
READ*, RFBL 

IF (RFBL.GT.@.@) GOTO 1780 
RFBL = 1.0D0 

WRITE (1,*) RFBL 

IF (HH.NE.6) GOTO 1900 
WRITE (6,181) 
+ ‘LENGTH (CM).'/? 
READ*, TFTL 

IF (TFTL.GT.@.%) GOTO 1830 
TFTL = 1.0D0 

WRITE (1,*) TFTL 

WRITE (6,186) 
a /? (DEFAULT = 
READ*, TFBL 

IF (TFBL.GT.@.0) GOTO 1880 
TFBL = 1.0 

WRITE (1,*) TFBL 

IF (HH.NE.7) GOTO 1990 
WRITE (6,191@) 

+ (DEFAULT = 
READ*, TRIBL 

IF (TRIBL.GT.@.8) GOTO 1930 


TRIBL = 1.0D0 


ee 


) 


) 


) 


) 


1.0)’ 


) 


WRITE (1,*) TRIBE 
WRITE (6,1941) 
FORMAT(/'TEST LINE THREE’ 


ENDEILE i 


WRITE(6,1991) © 
FORMAT(/'TEST LINE FOUR’ 
RETURN 

END 


) 


) 
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APPENDIX B 


PROGRAM FOR PLOTTING THE RADIATION PATTERN 


He Fe He HE He He FE HE HE HE FE TE Fe HE HE FE FE HE HE FE FE HE HE HE OE FE HE HE EE FE HE FE HE FE HE HE HE EE HHH HHH HHHHHHHHH HHH KH HH 


THIS PROGRAM IS TO PLOT THE POWER PER UNIT SOLID ANGLE RADIATED 
AT THE FREQUENCY NU. THE THEORY AND ALL EQUATIONS ARE DESCRIBED 
IN REF.3. IT ASSUMES THETA IN COL. 1 (DEGREES) AND POWER IN COL. 2 * 
(WATT) RESPECTIVELY. IT SUPPORTS A SINGLE PLOT OF MULTIPLE CURVES ~* 
AND USES THE COMPRS/DISSPOP OPTION OF DISSPLA. = 


Cw i ee A 


C 


¥* 


* * * * * 


REAL YX,XMIN,XMAX,YMIN, YMAX,X,Y,DX,AX,AY,XPHY,YPHY 

PNEBGER SLOUT , tamek, LTICKX, ITECKY, LYVAR, ICT, INN,NB,NETA,NF, 
+ NBL,LF,NMOD,NRAM,TLL,BL,NL 

CHARACTER*1 YES,ANS 

CHARACTER*4 BLNK(2) 

CHARACTER*6@ IMESS,IPHI,IRAM,IMOD,INB,IBL,IETA,LFL,INF,ITLL, 


+ INL 

C 

DIMENSION X(2000),Y(2000) 

DIMENSION YX(1:5000,1:2) 
C 

DATA YeS/’Y’/ 

DATA BLNK/’ eed cS 

Beo=.001 

C 
C INETIALISE PLOTTER 
C 


180 CONTINUE 
9400 FORMAT(A1) 
9401 FORMAT(1I1) 
9402 FORMAT(I2) 
WRITE(4,94@3 ) 
9403 FORMAT(' DO YOU WANT:’,/, 


X PeOOMP Rigs: yy | 
X EPENCOne: gu | 
X * ENTER NUMBER (11)’,/) 


READ(4,*) IOUT 

Pour EO.) THEN 
CALL COMPRS 

ELSEIF(IOUT.EQ.2) THEN 
CALL TEK618 

ENDIF 

AX = 4.6810 

AY = 6.7416 


T7 


CALL HWROT( ‘MOVIE’ ) 
CALL PAGE(8.5,11.@) 
XPHY=2.0840 
YPHY=2.7953 


DEFINE: LOWER/UPPER CASE STD AND GREEK SHIFT LABELS (+,!,?7,%), 
SHIFT FOR INSTRUCTION SET(&), AND THE TYPE FACE. 


MQ a0 4 


CALL MX1ALF(’STANDARD’,’!’) 

CALL MX2ALF(’L/CSTD’,’+’) 

CALL MX3ALF(’GREEK’,'??) 

CALL MX4ALF(’L/CGREEK’ ,’%’ ) 

CALL MX5ALF('INSTRU’,’&') 

CALL SWISSL 

CALL SHDCHR(90.,1,0.002,1) 
9410 FORMAT(D16.8) 

CALL HEIGHT(0.14) 

CALL PHYSOR( XPHY, YPHV) 

CALL AREA2D(AX,AY) 


C 
C CLOSE GRAPH AND REOPEN FOR LABEL AXES NUMBER TYPES 
C 
CALL ENDGR(@) 
CALL PHYSOR( XPHY, YPHY) 
CALL AREA2D(AX,AY) 
CALL XNONUM 
C 
C DETERMINE THE AXIS LABELS: 
C 
C X-AXIS LABEL AND NUMBER TYPE 
C 


CALL EXCMS('CLRSCRNV’ ) 
WRITE(4,9415) 
9415 FORMAT(’ DO YOU WANT INTEGERS ON THE X-AXIS? (Y/N)’,/) 
READ(4,9400) ANS 
IF(ANS.EQ.YES) CALL XINTAX 
CALL XNAME('A+NGLE! (?Q!, +DEGREES! )$’,10@) 


C Y-AXIS LABEL AND NUMBER TYPE 


WRITE(4,9417) 
9417 FORMAT(’ CHOOSE THE Y VARIABLE TO BE PLOTTED (11):’,/, 
x os W (WATT) ',/, 
X po ie oe DJ ua) 
READ(4,*) IYVAR 
IF (IYVAR.EQ.2) THEN 


CALL YNAME( ’D+J '$’,100) 
ELSE 

CALL YNAME('’W (W+ATT! )$’ ,180) 
ENDIF 


C 
C ROTATE NUMBERS ON Y-AXIS AND REVERSE TICKS 


te 


©2 


2 


aqQaqaaq 


790@ 


9425 


94.24 


9425 


poe ( 


9426 


9429 


CALL YAXANG(@. ) 

CALL XREVTK 

CALL YREVTK 

WRITE(4,790@) 

BORA sNeUL YMIN OF THE FIRST SUBPLOT (F16.8):’,/) 
READ(4,*) YMIN 

CALL YAXEND ('NOLAST’ ) 
YMAX=YMIN+.1 

Probe CRAN CO. 1.4 1.,YMIN, YMAX,. 1): 
CALL THKFRM( .927) 

CALL FRAME 

CALL ENDGR(@) 


CLOSED GRAF FOR AXIS LABELS, NOW REOPEN TO PLOT 


CALL PHYSOR( XPHY, YPHY ) 
CALL AREA2D(AX,AY) 
CALL RESET ('YAXEND' ) 
CALL RESET ('XNONUM’ ) 
CALL YNAME( BLNK,2) 
CALL XNAME(BLNK,2) 
CALL ENDGR(@) 


POSITION OR REPOSITION ORIGIN AND SCALE PLOT 


CALL PHYSOR( XPHY, YPHY) 
CALL AREA2D(AX,AY) 
CALL RESET ('YAXEND’ ) 
CALL YAXEND('NOFIRST’ ) 


DEFINE COORDINATE SYSTEM 
KEEP OR CHANGE CONSTANTS 


CALL EXCMS('CLRSCRN’ ) 
WRITE(4,9423) 

FORMAT(’ INPUT XMIN: (F16.8)’) 
READ(4,*) XMIN 

WRITE(4,9424) 

FORMAT(’ INPUT XSTEP: (F16.8)’) 
READ(4,*) XSTEP 

WRITE(4,9425 ) 

FORMAT(’ INPUT XMAX: (F16.8)') 
READ(4,*) XMAX 

WRITE(4,9427) 

FORMAT(’ INPUT YMIN: (F16.8)’) 
READ(4,*) YMIN 

WRITE(4,9426) 

FORMAT(’ INPUT YSTEP: (F16.8)’) 
READ(4,*) YSTEP 

WRITE(4,9429) 

FORMAT(’ INPUT YMAX: (F16.8)’) 
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©2 


9330 


See 


READ(4,*) YMAX 

CALL EXCMS('CLRSCRN’ ) 

WRITE(4,9330) 

FORMAT(’ HOW MANY TICKS PER X-AXIS STEP? (I2)’) 
READ(4,*) ITICKX 

CALL XTICKS(ITICKX) 

WRITE(4,9331) 

FORMAT(’ HOW MANY TICKS PER Y-AXIS STEP? (I2)’) 
READ(4,*) ITICKY 

CALL YTICKS( ITICKY ) 

IF(XSTEP.EQ.0.0) XSTEP=(XMAX-XMIN)/AX 
IF(YSTEP.EQ.0.0) YSTEP=(YMAX-YMIN)/AY 

CALL GRAF(XMIN,XSTEP,XMAX,YMIN, YSTEP , YMAX ) 

CALL YNONUM 

CALL YGRAXS(YMIN,YSTEP,YMAX,AY,’ $',-100,AX,0.@) 
CALL RESET ('YNONUM' ) 

CALL XNONUM 

CALL XGRAXS(XMIN,XSTEP,XMAX,AX,’ $',-100,0.0,AY) 
CALL RESET ('XNONUM'’ ) 


CURVE PLOTTING 
CALL LINESP(1.8) 
CALL SPLINE 
CALL THKCRV(0.02) 
CALL GRACE (0.980) 


DRAW A CONNECTED CURVE 


MULTIPLE CURVE RETURN POINT 


2000 


9445 


9455 


CONTINUE 
CALL EXCMS('CLRSCRN'’ ) 
WRITE(4,9445) 


FORMAT(’ YOU HAVE A CHOICE OF LINE TYPE:’,/, 

+ , {eos Soc eSeae SOLID LINE’',/, 

+ ; ir mre et ee es DOTTED LINE’ ,7; 

i ’ ee ee DASHED LINE’,/, 

+ fi cle a eens CHAIN DOT LINE’,/, 

in ee eee ee CHAIN DASH LINE’,/, 

+ ' PLEASE INPUT LINE TYPE NUMBER: (11)',/) 


READ(4,*) ICT 
IF(ICT.EQ.2) CALL DOT 
IF(ICT.EQ.3) CALL DASH 
IF(ICT.EQ.4) CALL CHNDOT 
IF CLOT EO Sie caAbb CHNDSH 
CALL EXCMS('CLRSCRN' ) 


DATA INPUT FROM EXTERNAL DATA FILE 
WRITE(4,9455) 


FORMAT(’' INPUT THE FILE DATA SET NUMBER (I2)') 
READ(4,*) INN 


80 


9460 


2400 
8000 


8019 


Q 


Q 


9475 


ome i? i®) 


9480 


94835 


NUCOLS = 2 
IN=@ 
DO 2400 II=1,5000 
READ(INN,9460) (YX(II,JJ),JJ=1,NUCOLS) 


READ (4,9401) NF 

IF (NF.EQ.1) THEN 
INF= ’GAUSSIAN FUNCTIONS’ 
WRITE (4,9484) 


81 


FORMAT(4X,F8.3,2X,F20.18) 
thee. ied) bis] O))-ANDe( [1.GT.2)) GO TO 8000 
IN=IN+1 

CONTINUE 

CONTINUE 

ICOL = 2 

DO 8019 II=1,I1N 

Mell =y¥ x LE. 1°) 

Yl J=YxXC EE, ECOL) 

CONTINUE 

REWIND INN 
TERMINATE DATA READING & PLOT CURVE 

CALL CURVE( X,Y,IN,9) 
RESET LINE TYPE 

ieClCL.s60-2 ) CALL RESET( *DOT’ ) 

Le Crentsne, >) CALE RESET *DASH’ ) 

IF(ICT.EQ.4) CALL RESET( ’'CHNDOT’? ) 

IF(ICT.EQ.5) CALL RESET( 'CHNDSH’ ) 

CALL EXCMS( 'CLRSCRN’ ) 

WRITE(4,9473) 

FORMAT(’ DO YOU WISH TO TERMINATE THIS SUBPLOT?’ ) 

READ(4,9400) ANS 

IF(ANS.NE.YES) GO TO 2000 

END OF PLOT 

WRITE ADDITIONAL LABELS 

CALL ENDGR(@) 

CALL AREA2D(AX,AY) 

CALL EXCMS( 'CLRSCRN’ ) 

WRITE(4,9480) 

FORMAT(’ DO YOU WANT ADDITIONAL LABELS ?’,/) 

READ(4,9400) ANS 

IF(ANS.EQ.’N’) GO TO 9900 

WRITE(4,9483) 

FORMAT(’ YOU HAVE A CHOICE OF FUNCTION :',/, 

+ ; (a GAUSSIAN FUNCTION’,/, 

+ ———— LEVEL FUNCTION’, /, 

+ o==—---=— TRAPEZOIDAL FUNCTION’ ,/, 

+ * INPUT THE NUMBER OF FUNCTION TYPE (11) :’,/) 


9484 


9486 


+++ +++ + + + 


+++ ++ 4+ + + + 


FORMAT( ’SPECIFY BUNCH LENGTH (F4.2):',/, 


== 0.24 cM ',/, 
, = 0.88 CM ’,/, 
’ 3 S22S5 Treo CM: {aye 
’ / ee Se 1.76 Chee 
5 sas 2.40 7CNe ee = 
= 2.63 CM. aes 
, ———— 3.50. ‘Chime 
, St22--2 A038. CMs 
‘INPUT THE NUMBER(I1):___',/) 


READ(4,9401) NBL 
IF (NBL.EQ.1) THEN 
IBL= ’BL = 0.24 +CM!$’ 
ELSEIF(NBL.EQ.2) THEN 
IBL= ‘BL = 90.88 +CM!$’ 
ELSEIF(NBL.EQ.3) THEN 
IBL= ‘BL = 1.20 +CM!$’ 
ELSEIF(NBL.EQ.4) THEN 
IBL= ‘BL = 1.76 +CM!$’ 
ELSEIF(NBL.EQ.5) THEN 
IBL= *’BL = 2.40 +CM!$’ 
ELSEIF(NBL.EQ.6) THEN 
TBieer Ble) 2.60587: SD - 
ELSEIF(NBL.EQ.7) THEN 
TBL=: *BEw= "3.90 4CM' 3" 
ELSE 
IBL= *BL = 4.38 +CM!$’ 
ENDIF 
ELSEIF(NF.EQ.2) THEN 
INF= 'LEVEL FUNCTIONS’ 
WRITE(4,9486) 
FORMAT( ’SPECIFY LEVEL FUNCTI 


ON 


LENGTH (F422): jae 


’ jones 1200 =CMiauia 
' a 5.00 Ch ae 
, Sj Sao AA SOmON ee 
: 4 ----- 5 O0n CMa ay. 
5 Sees = 6200 CM way. 
; 6 700 gCMe ee 
’ ee 8.008 CM ',/, 
, Ss fue 9700 CM * 7; 
‘INPUT THE NUMBER( 11): aay) 


READ(4,9401) LF 
IF oCLE. EO. 1) THEN 

LFL= ’LFL = 1.800 +CM!$’ 
ELSEIF(LF.EQ.2) THEN 

LFL= 'LFL = 35.00 +CM!$’ 
ELSEIF(LF.EQ.3) THEN 

LFL= 'LFL = 4.50 +CM!$’ 
ELSEIF(LF.EQ.4) THEN 

LFL= *LFL = 5.00 +CM!$’ 
ELSEIF(LF.EQ.5) THEN 

LFL= 'LFL = 6.00 +CM!$' 
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9501 


++ ete ttt t 


9502 


+++ + 4+ 


ELSEIF(LF.EQ.6) THEN 
LFL= "LFL = 7.99 +CM!3’ 
ELSEIF(LF.EQ.7) THEN 
LFL= 'LFL = 8.99 +CM!$' 
ELSE 
LFL= 'LFL = 9.90 +CM!$’ 
ENDIF 
ELSEIF(NF.EQ.6) THEN 
INF= 'TRAPEZOIDAL FUNCTIONS’ 
WRITE(4,9591 ) 


FORMAT('SPECIFY TOP LEVEL LENGTH 

1 ---- 9.97 CM 
-94 CM 
plea cl! 


.2 CM 
9 CM 


TIA WT FP UW fo 
ou FW fh 


READ(4,9401) TLL 
IF (TLL.EQ.1) THEN 
ITLL= *TFTL = 8.07 +CM!$’ 
ELSEIF (TLL.EQ.2) THEN 
TPEL= (TFTL = 2.94 +CM!$' 
ELSEIF (TLL.EQ.3) THEN 
PEG TPETL = 35.185. +CM!S$' 
ELSEIF(TLL.EQ.4) THEN 
TTob= “TTR = 420° 4CM!$’ 
ELSEIF(TLL.EQ.5) THEN 
ITLL= ‘TFTL = 5.9 +CM!$’ 
ELSEIF(TLL.EQ.6) THEN 
ITLL= ‘TFTL = 5.04 +CM!$’ 
ELSE 
ffbe- TERETE = 8.96 +CM!$' 
ENDIF 
WRITE(4,95@2) 
FORMAT(’SPECIFY BASE LENGTH :' 


-94 CM 
-96 CM 
’ INPUT THE NUMBER (11): uae) 


ey ae 
ae ae 
om 
a. 

An ike 

Se an 
4/, 


sls 


; Ie ==-= 5.215 CM” ,/. 


2 ---- 7.98 CM 
3 ---- 8.9 CM 
; 4 ---- 8.96 CM 


' INPUT THE NUMBER (11): re ie, 


READ(4,9491) BL 
Pe (Bien. 1) fHEN 

PBib=) “TFBL = 5.215 +€M!3$’ 
ELSEIF(BL.EQ.2) THEN 

IBL= 'TFBL = 7.9 +CM!$’ 
ELSEIF(BL.EQ.3) THEN 

IBL= ‘*TFBL = 8.9 +CM!$’ 
ELSE 

IBL= ‘'TFTL = 8.96 +CM!$’ 
ENDIF 
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a ae 


ae 


oda ae 


ENDIF 


WRITE( 4, 9485 ) 


9485 FORMAT(' 


SPECIFY RADIATION MODE (Ti 7] 


su 1===— CERENKOV = 
+ ’ 2 --- SUB-CERENKOV'’',/, 
+ ' INPUT THE NUMBER: __',/) 
READ(4,9401) NMOD 
IF(NMOD.EQ.1) THEN 
IMOD = °'C+ERENKOV!S! 
ELSE 
IMOD = '+SUB-!C+ERENKOV!$' 
ENDIF 
WRITE(4,9487) 
9487 FORMAT(' SPECIFY ENERGY OR NBETA (I1):’,//, 
X i FORGE (476 
X 1 -- E = 2.555 MEV’ ,/, 
X ' FOR NBETA’,/, 
X 2 -.-- NBETA = @.99 NM’,/, 
X , 3 -- NBETA = 0.98 NM’,/, 
X 4 -- NBETA = 0.97 NM’,/, 
X , 5 -- NBETA = 0.96 NM’,/, 
X ’ 6 -- NBETA = 0.95 NM’,/, 
X ' INPUT THE NUMBER: _',/) 
READ(4,9401) NB 
IF(NB.EQ.1) THEN 
e INB = '+N!%B! = 1.0098’ 
INB = 'E = 2.555 M+E!V3S’ 
ELSEIF(NB.EQ.2) THEN 
INB = '+N!%@Bt = 0.998’ 
ELSEIF(NB.EQ.3) THEN 
INB = °+N!%Bt = 0.988” 
ELSEIF(NB.EQ.4) THEN 
INB = °+N!@B! = 0.978’ 
ELSEIF(NB.EQ.5) THEN 
INB = ’+N!%@B! = 0.968’ 
ELSE 
INB = '+N!%B!t = 0.958’ 
ENDIF 
CALL EXCMS('CLRSCRN’ ) 
WRITE(4,9488) 
9488 FORMAT(’ SPECIFY ETA (I1):’,//, 
X : | =="BTA = 1257 37- 
X Oo SORT N eso 
x 3 =e BE os ewes 
X ' 4 -= BTA = 74 32 
X ’ INPUT THE NUMBER: __',/) 
READ(4,9401) NETA 
IF(NETA.EQ.1) THEN 
IETA = '4CP =912 $) 
ELSEIF(NETA.EQ.2) THEN 
IETA = '%C! = 24 $° 


ELSEIF(NETA.EQ.3) THEN 
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9490 


PS PS PS OOS PS ORS 


PEA = eC! 
ELSE 

IETA = ‘SC! 
ENDIF 


WRITE(4,9490) 
SPECIFY BEAM PATH LENGTH (I11):’,//, 


FORMAT ( ' 


=r aD 


= 74 $’ 


1 == L = 209CM (BETA=12,N=6)’,/, 
2 -- LZ = 42 CM (ETA=12,N=3)’,/, 
pa = e648 CM (ETA=37,N=6)",/, 
Ave) 126.2 CM (ETA=12,N=1)",/, 
; pe oem Oe ome ( ETA=57,N=5)°,/, 
, INPUT THE NUMBER: _ ',/) 

READ(4,9401) NL 

IF(NL.EQ.1) THEN 

ieieee—nt y= 2] CM tS 
ELSEIF(NL.EQ.2) THEN 

INL = 'L = 42 +CM!$’ 
ELSEIF(NL.EQ.3) ‘HEN 

INL = 'L = 64.8 +CM!$’ 
ELSEIF(NL.EQ.4) THEN 

INL = ’L = 126.2 +CM!$’ 
ELSE 

INL = 'L = 129.5 +CM!$’ 

ENDIF 

WRITE(4,9492) 

FORMAT(’ SPECIFY WAVE LENGTH (I1):’,//, 

1 -- RAMDA = 10.52 CM (N=1)’,/, 
2 -- RAMDA = 3.50 CM (N=3)’,/, 
5 e== RAMDAS= 1-7 50GM (N=6)',/, 
’ INPUT THE NUMBER: _ ',/) 
READ(4,9401) NRAM 
IF(NRAM.EQ.1) THEN 
IRAM = ‘SL! = 10.52 +CM!$’ 
ELSEIF(NRAM.EQ.2) THEN 
PeAMio= "4! = 5.50 +CM'S’ 
ELSE 
IRAM = '%L! = 1.75 +CM!$? 

ENDIF 

XMES = 2.880 

YMES = 6.130 

IMESS = IMOD 


CALL MESSAG( IMESS, 100,XMES, YMES) 


YMES = 


IMESS = INB 


YMES - 0.2590 


CALL MESSAG( IMESS,100,XMES, YMES) 


YMES = 
IMESS = 


YMES - @.25@ 
IETA 


CALL MESSAG( IMESS,10@,XMES, YMES) 
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9493 


9495 


9999 


YMES = YMES - 9.2590 
IMESS = INL 
CALL MESSAG( IMESS,19@,XMES, YMES ) 


YMES = YMES - @.2590 
IMESS = IRAM 
CALL MESSAG( IMESS,190,XMES, YMES) 


YMES = YMES - 9.259 
IF(NF.EQ.1) THEN 
IMESS=IBL 
ELSEIF(NF.EQ.2) THEN 
IMESS=LFL 
ELSEIF(NF.EQ.6) THEN 
IMESS= ITLL 
CALL MESSAG( IMESS,100@,XMES, YMES) 
YMES = YMES-@.259 
IMESS = IBL 
ENDIF 
CALL MESSAG( IMESS,19@,XMES, YMES ) 


IF(HH.EQ.6) THEN 
XMES = 1.29 
ELSE 
XMES = 1.50 
ENDIF 
YMES = -2.25 
IMESS = INF 
CALL MESSAG( IMESS,199,XMES,YMES) 


CONTINUE 
CALL EXCMS('CLRSCRN’ ) 
IF(IOUT.EQ.2) THEN 
WRITE(4,9493) 
FORMAT(' <CLEAR>, THEN READY FOR HARD COPY’ ) 
ENDIF 
CALL ENDPL(0) 
WRITE(4,9495) 
FORMAT(' DO YOU WANT TO MAKE ANOTHER FIGURE?’,/) 
READ(4,9400) ANS 
IF(ANS.NE.YES) THEN 
GO TO 9999 
ELSE 
ENDIF 
CALL RESET(’ALL’ ) 
GO TO 100 


CONTINUE 
CALL DONEPL 
Sy bl 6) & 

END 
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